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ON THE INTEGRAL KERNELS OF DERIVATIVES OF THE 
ORNSTEIN-UHLENBECK SEMIGROUP 

JONAS TEUWEN 


Abstract. This paper presents a closed-form expression for the integral ker¬ 
nels associated with the derivatives of the Ornstein-Uhlenbeck semigroup e^^. 
Our approach is to expand the Mehler kernel into Hermite polynomials and 
applying the powers of the Ornstein-Uhlenbeck operator to it, where we 
exploit the fact that the Hermite polynomials are eigenfunctions for L. As an 
application we give an alternative proof of the kernel estimates by m, making 
all relevant quantities explicit. 


1. Introduction 


Much effort [U [31 IH 0 El 13 HI HHl (HI (HI 113] has gone into developing the 
harmonic analysis of the Ornstein- Uhlenbeck operator 

(1) L:=^A-{x,V). 

On the space L'^ , d'y), where 7 is the Gaussian measure 

(2) d7(a:) := dx, 


this operator can be viewed as the Gaussian counterpart of the Laplace operator 
A. Indeed, one has L = — V*V, where V is the the usual gradient and the V* is its 
adjoint in L^(R'^, d 7 ). It is a classical fact that the semigroup operators t > 0, 
are integral operators of the form 

e*^u(-) = / Mt{-,y)u{y)d'y{y), 
where Mt is the so-called Mehler kernel [12] 


( 3 ) 


Mt{x,y) = 


exp 


|e ‘x-yl" 

-2t 


1 — e 


(1 _e-2t)d/2 




(see [13] for a representation of Mt which makes the symmetry in x and y explicit). 

Developing a Hardy space theory for L is the subject of active current research 
iniis]- In this theory the derivatives {d^/dt^)e*^ = L^e*^ play an important role. 
The aim of the present paper is to derive closed form expressions for the integral 
kernels of these derivatives, that is, to determine explicitly the kernels such 
that we have the identity 


( 4 ) 


lW)=/ Mt^{;y)u{y)d^{y). 


Direct application of the derivatives d^/dt^ to the Mehler kernel yields expressions 
which become intractible even for small values of N. Our approach will be to 
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expand the Mehler kernel in terms of the L^-normalised Hermite polynomials and 
then to apply to it, thus exploiting the fact that the Hermite polynomials are 
eigenfunctions for L. 

As an application of our main result, which is proved in section [Rafter developing 
some preliminary material in the sections [21141 we shall give a direct proof for the 
kernel bounds of m in section n 


2. Preliminaries 


In this preliminary section we collect some standard properties of Hermite poly¬ 
nomials and their relationship with the Ornstein-Uhlenbeck operator. Most of this 
material is classical and can be found in [Hill]. 


2.1. Hermite polynomials. The Hermite polynomials u > 0, are defined by 
Rodrigues’s formula 

(5) i?„(x):=(-ire"'a:e-"'. 

Their L^-normalizations, 


hrj. ■ — 


Hr, 


V2’^n! 


form an orthonormal basis for T^(R, dy). We shall use the fact that the generating 
function for the Hermite polynomials is given by 


( 6 ) 


E Hnix) ^ri _ 2tx-t^ 


n—0 


The relationship with the Ornstein-Uhlenbeck operator is encoded in the eigenvalue 
identity LHn = —nH^, from which it follows that for all t > 0 we have = 

e“*"R„. From this one quickly deduces that the Mehler kernel is given by 


(7) 




n—0 


We will need two further identities for the Hermite polynomials which can be found, 
e.g., in [TJ Chapter 18]: the integral representation 


( 8 ) Hnix) = 

and the “binomial” identity 

(9) 


{-2i)^e^ 




Ce 




Hnix + y) = ^(fj{2yr->^Hkix). 

k=0 ^ ' 


2.2. Hermite polynomials in several variables. The Hermite polynomials on 
are defined, for multiindices a = (oi,..., aa) G by the tensor extensions 

d 

( 10 ) Ha{x) := hatiXi), 

n—1 

for X = {xi ,..., Xd) in R'^. The normalized Hermite polynomials 

(11) ha := , where a! = Oi!. 

\/2l“la! 
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form an orthonormal basis in L^(R'^,d 7 ), and we have the eigenvalue identity 

(12) LHa = —\a\Ha, where |a| = oi + • • • + ad- 

If we consider the action of L^e*^ on a Hermite polynomial through the multi¬ 
nomial theorem applied to \a\^ we get (writing Ld for the operator L in dimension 
d and Li for the operator L in dimension 1) 

(13) 

L^e*^ha{x) = |a|^e (ici) •... • ha^{xd) 

= XI f ^ • ... • •... • e“‘^/iQi(a;i) • ... •/iad(a;d) 

|„|=Ar \ni,n2,...,ndj 

= X f ^ Xi^e‘^^/iai(a;i)-...-Li‘^e‘^^/ia^(a;d)- 

This implies that we can reduce the question of computing the d-dimensional version 
of the integral kernel to the one-dimension one. 

3. A combinatorial lemma 

From now on we concentrate on the Ornstein-Uhlenbeck operator L in one di¬ 
mension, i.e., in (R,d7). We are going to follow the approach of [T^. Recalling 
the identity L/i„ = —nhn, we will apply to the generating function of the Her¬ 
mite polynomials ([5]). A problem which immediately occurs is that A and (x, V) 
do not commute, and because of this we cannot use a standard binomial formula 
to evaluate . Instead, we note that 

Lg = -tdtg. 

In particular this implies that 

L^g = i-l)'^DNg, 

where 

(14) Dn ■■= tdt otdto ■■ ■ o tdt_ = 

k times 

The following lemma will be very useful. 

Lemma 1. IFe have 

N f ^ 

( 15 ) 

n=0 ^ ^ 

where are the Stirling numbers of the second kind. 

The Stirling numbers of the second kind are quite well-known combinatorial 
objects. For the sake of completeness we will state their definition and recall some 
relevant properties below. For more information we refer the reader to |15j . The 
related Stirling numbers of the first kind will not be needed here. 

We begin by recalling the definition of falling factorial 

(16) (j)n := j{j -l)...{j-n + l)= , 

(j - n)'. 


for non-negative integers k > n. 
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Definition 1. For non-negative integers N > n, the number Stirling number of 
the second kind {^} is defined as the number of partitions of an N-set into n 
non-empty subsets. 


These numbers satisfy the recursion identity 


(17) 



= N 


N-l 

n 


N-l 
n — 1 


For all non-negative integers j and k one has the generating function identity 

N ^ y. 

( 18 ) = 

n=Q '' t 


4. Weyl Polynomials 

Before turning to the proof of lemma [U let us already mention that it only 
depends on the commutator identity [t, dt\ = —1. This brings us to the observation 
that Weyl polynomials provide the natural habitat for our expressions. Rougly 
speaking, a Weyl polynomial is a polynomial in two non-commuting variables x 
and y which satisfy the commutator identity {x,y\ = —1. This is made more 
precise in the following definition. 

Definition 2. The Weyl algebra over a field F of characteristic zero is the ring 
F(x, y) of all polynomials of the form p{x, y) = X)m=o Z)n=o ^mnx'^y'^ with coeffi¬ 
cients Cmn G F in two noncommuting variables x and y which satisfy the commu¬ 
tator identity 

[x,y] := xy-yx = -1. 

We now have the following abstract version of lemma [T] 

Lemma 2. In the Weyl algebra F{x,y) we have the identity 

™ r 1 

( 19 ) (xy)"' = |7 

where {™} are the Stirling numbers of the second kind. 

As a preparation for the proof of lemma[2]we make a couple of easy computations. 
If we set D := xy, then 

( 20 ) Dx^ =x^D + mx^, 

( 21 ) Dy^ = y”^D - my^. 

This can be shown by induction on m. For instance, note that 
Dx™ = x{D + l)x™-i = xDx^-^ + x’". 

If we take this a bit further and have p G F[Z?], then 

(22) p{D)x"^ = x"^p{D+ m), 

(23) piD)y^ = y^p{D-m). 
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The TO-th powers, m > 1, of a; and y satisfy 


(24) 

m—1 

= n (^ - *), 
2=0 

m 

(25) 

y'^x^ = \{{D + i). 
2=1 


This can be seen using induction: 

and 

^m+l^m+l ^ Ip + 1)), 

Definition 3. The weighted degree of a monomial a;™y” G F{x,y) is the integer 
m — n. A polynomial in F{x^y) is said to he homogeneous of weighted degree j if 
all its constituting monomials have weighted degree j. 

Left multiplication by xy is homogeneity preserving, i.e., for all j G Z it maps 
the set of homogeneous monomials of weighted degree j into itself. To prove this, 
first consider a monomial x'^y'^ of dweighted egree j = m — n. Then, 

{xy)x"^y^ ^ {x'^[xy) + mx'^)y'^ = x'^{xy)y^ + mx^y^ = x’^+^y'^+^ + mx^y^, 

and we see that weighted degree of homogeneity is indeed preserved. The general 
case follows immediately. Through (I24|) we conclude that left multiplication 
is homogeneity preserving as well, for all non-negative integers k. We claim that 
left multiplication by x'^y^ is homogeneity preserving only if i = j. To see this note 
that 

yx'^y^ = x^y^~^^ + ix^~^y^ 
from which we can deduce that 

x'^y^x'^y^ = -h lower order terms. 

From which the claim follows. 

Finally, a polynomial is homogeneity preserving if and only if all of its constitut¬ 
ing monomials have this property. If this were not to be the case we could look at 
the highest-order non-homogeneous term and note from above x'^y^x^y^ would 
give terms of a lower order in the polynomial expansion which cannot cancel as 
they have different powers of a; or y. 

It follows from these observations that 

r ^ 

(26) ^0 := { E 

^ n=0 

is precisely the set of homogeneity preserving polynomials in F{x,y). 

Now everything is in place to give the proof of lemma [51 

Proof of lemma[^ As (xy)^ is homogeneity preserving, we infer that there are co¬ 
efficients in F such that 

k 

ixy)^ = J2a^xY. 

2=0 



( 27 ) 
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We will apply to the right on both sides of (j27|l and derive an expression for the 
af. First note that (l22l) gives 

{xyfx^ = x^{xy + jf, 
and ((24)l together with ((2^ gives 

_ i—1 _ j—1 

x''y''x^ = — i)x^ = x^ ]^(a:y — i + j). 

i=o e=o 

Hence, to find the coefficients of it is sufficient to consider 

k i—1 

{xy + jf = '^al^nW{xy-(. + j). 

i=0 1=0 

Comparing the constant terms on both the left-hand side and right-hand side, we 
find 

k 2—1 k 

( 28 ) = = 

2—0 i—0 2—0 

where {j)i is the falling factorial as in (fTCl) . Comparing (1^51) with the generating 
function of the Stirling numbers of the second kind {f} as given in (fTSl) . we see 
that of = {f}. This concludes the proof of lemma [5] □ 

5. The integral kernel of 

As mentioned before, as a first step we would like to apply to the generating 
function g{x,t) := for the Hermite polynomials ([n]). We 

first compute the action of on the generating function. 

Lemma 3. We have 

(29) QN^-(x-tfW ^ Hn{x - t). 

Proof. We first note that, 

= 2{x - t)e-(^-*)' = 

Using this we get 

= (_)2g.^gN-2^2g-(.-i) = 


= (-l)^e“'(9fe-(^-*)'. 

By a change of variables, 

V=x-t 

Hence, by (O, 

qN^-(x-W+x^ = - t). 

□ 
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Lemma 4. For all x G R and t > 0 we have 


(30) 


n=0 


Proof. This is now easy to prove. Recalling that L = —tdt and using (fT51) . we get 


= (-1)^ 

n=0 


N 


= (-l)^e-("-‘) ^ I ^ - t). 


n—0 


N 


□ 


Our next theorem is the main result of this paper and provides an explicit ex¬ 
pression for the integral kernel of 

Theorem 1. Let L be the Ornstein-Uhlenbeck operator on L^(R‘^,d 7 ), let t > 0, 
and let N > 0 be an integer. The integral kernel of L^e*^ is given by 


M^ix,y) = Mtix^y) 

(31) l"l=^ 


N 

Til J ■ ■ ■ 5 


rii rrii 


nEE2 

i=l mi=0 £i=0 


Hi I /mi 


mi 


X - 


\/l — e 


2t 


Irrii-li 


Hi- {Xt)H2mi-ii 


XjB * - Pi 

\/l — e“^‘ 


Proof. We first prove the result for d = 1. Pulling through the integral expres¬ 
sion o for involving the Mehler kernel, we must find a suitable expression for 
the kernel Mt{-,y). Using ([7]) and the normalization of Hm in (ED 

we get 


M,^ix,y)=L^Y 

^ 'Y 




•m 


ml 2 
e”*"'" 1 


Hjn (jj) 


m—0 

ry—tm 


m 


= L 


11 


N ' 


m—0 


! 2 


■m 


i^„,(x)^^e 


m ^ poo 


e-« dC 


e-« 


ml 


m—0 


N ' 




df. 


The operator is applied with respect to x here, so by lemma 0] we get 

y'^ pOO 

MYx,y) = ^ / d^ 

V7-00 
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where in last line we have used the analytic continuation of the algebraic identity 
(|30ll . Similarly we can expand Hm{y + using ([9|). This gives 


so that MF can be written as 
(32) 


X -\-y ( at') / \ /*oo 

If we set M = 2m — i, this reduces our task to computing the integral 




(33) 


(*ee-*)" de 

pOO 

V ^ J — oo 


To make the computation less convolved, let us set 

/ - ^ — y 

at-- VI- e-2*, and, y) := -j===. 

This allows us to write the exponential in the integral (13311 as 

This reduces the problem, after the substitution at^ —> to computing the integral 


aV — /-OO 

® ® ' Q2il3tix,y)^^-^-‘ ^ 


rz: „M+i 
V ^ at j —oo 

The final integral is an integral expression for the Hermite polynomials ([5]) , so 


.oo 










Next we note that exp(j/^ — (3t{x, yY)at ^ = Mt, the Mehler kernel from (|3]). Hence, 


N 


(34) 


MVix,y) = Mt{x,y) ^ ^ 


m\ f N 


X Hl(x)H2m-l 


m —0 £—0 

( xe~* — y 


Cl — e 


2t 


2 m—£ 


V Cl — e“^* 

Applying (11^1) we get the result in d dimensions. 


□ 
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6. An application 

As an application of our main result, in this section we give an alternative proof 
of the bounds on the kernels K and K of m (see the definition below), making 
the dependence on the parameters more explicit. These kernels play a central role 
in the study of the Hardy space i?^(R‘^,d 7 ) in [11], where the standard Calderon 
reproducing formula is replaced by 

(35) u = C f (t^L)^^^et^M—+ [ rtdy, 

JQ t v/R-i 

where C is a suitable constant only depending on N and a (this can be seen by 
letting u be a Hermite polynomial). The kernels K and K then occur in several 
decompositions, and the estimates below allow them to be related to classical results 
about the Mehler kernel. 


Definition 4. We define the kernels K and K by 

[ y)u{y) d'y{y) = (f L)^^u{x), 

v/Rd 

f Kt^,N,a,jix, y)u{y) d 7 (y) = (t^L)^e^^td* u{x). 


Note that the operators on the right-hand sides are indeed given by integral 
kernels: the first is a scaled version of the operator we have already been studying, 
and a duality argument implies that the second is given by the integral kernel 

Kt2^N,a,j{x,y) = tdxjKt2^N^a{x,y). 

Thus, both kernels are given as appropriate derivatives of the Mehler kernel. 

We begin with a technical lemma which is a rephrased version of El Lemma 
3.4]. One should take note that we define the kernels with respect to the Gaussian 
measure whereas, El defines these with respect to the Lebesgue measure. 


Lemma 5. For all a > 1 and all t > 0 and x, y in we have 


"ex- yp 


(36) 

1 — e 

Additionally, we have 

(37) 


> Tre 


-2t 


|e —min (|xf, |yf) 


1 — e 


-2t 


1 — e 


ae < 


1 — e 


-2t 


1 — e 


< a. 


Theorem 2. Let N be a positive integer, 0 < t < T. The for all large enough 
a > 1 we have 

(1) Ift\x\ < C, then 

\Kt-2,a^N{x,y)\ <T,N aexp{^C'^)Mt2{x,y)exp(^--^^ )' 

(2) Ift\x\ < C, then 

\Kt2^a,NA^,y)\ <T,N aexp(^C'^)Mt2(a:,?/)exp^-^^i^j-— 
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Proof. For we use Theorem [T] to obtain, after taking absolute values, 

d rii rrii 

\Kt 2 ^a,N{.x,y)\< M^{x,y) ^ 

|fe|=iV 


e “ 


N 

Tl\ , . . . , TlfiJ ^ 
2mi-ii 




■mi I 

e,. 


\Hii {xi)\ 


i—1 ^i—Omi—Q 

-£ 

XiB “ - y^ 


m, \ [ n, 
m, 


H2mi—ti 


1 — e 


-2-^ 


'l-e-2^- 

Recalling that ii-\ - + id < N, using the assumptions t <T and t|a;| < C we can 

bound by considering the highest order term to obtain 

^2fc. I 

Using (1571) we proceed by looking at 

-2p \ 1/2 — 2/p' 


M^{x,y) = Mt 2 {x,y) 


l-e" 


exp 


< aMt 2 (x, y) exp 


1 - e-2^ 

e-‘'cr-yp 


1 - e-2‘= 


exp 


e “ X 


-2/P 


1 -e- 


1 - e-2*" 


exp - 


1 |e “X-i/p 


1 — e 


We can now bound the final Hermite polynomial in the expression of the kernel. 
Setting Mi = 2mi — £{ we get 


1 — e 


-2^ 


Mi 


H, 


Mi 


e ^Xi-yi 


1 — e 


-2^ 


< 


N 


-2^ 


< 


1 — e 


|e ^ x^- y^ 


Mi _t± \ \ Mi 

e ^Xi-yi\ 


1 — e 


-2-^ 


1 — e 


Mi 


Also, 


|e ^Xi-yi\ 

1 -2^ 

I - e ^ o. 


Mi 


exp - 


1 |e 

1 - e 


< 1 . 


Putting these estimates together, using Lemma [SJ and taking a > 1 so large that 


a a np. P 

1 —rw < — TTW 3'iid 1 — e “ > —. 
4e2^ “ “ a 


we obtain 


\Kf^,a.N{x,y)\ <T,N 0!Mt2{x,y)exp 


\e * X- 2/p 
1 - e-2*= 


exp - 


< aM (2 (x, y) exp 


|e * X- 2/P 


1 - e-2*= 


exp - 


1 | e °x- 2 /p 
1 - e 
a |e“*^x —2/P 
4e^^ 1 — e“^*^ 


exp 


1 C|xp 

2 1 -2-R 

^ 1 — e “ 


= aMt2(x,2/)exp^(^ 
< q;M (2 (x, 2 /) expl — 


1 - 


a N |e ‘ X — 2/p 


4e^^/ 1 — e“^* 
a |e“‘ X — 2/p\ (a 


1 C|xp 

57773? 


8e^^ 1 — e“2‘^ 


< Of exp ( —C^ )Mt 2 (x, 2 /) expl — 




|e * X - 2 /p 


8e^'^ 1 — e“2*^ 
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using t\x\ < C in the last step. 

For the bound on K we consider 


tda- 




XjB * - Hi 

\/l — e“^‘ 


= tH„ 


XjB * - Hi 

‘ V Vl — e“^‘ 


+ Hi.{xi)td^,Hr, 


dxi^liXi) 

( Xie~* - Ui 


VVl-e-2‘, 

So, as the first term on the right-hand side just decreases in degree we look at 


td:r. 


XjB * - Vi 

\/l — e“^‘ 


\ mi —1 

XiB - y, 


y/l — e“^ 


\/l — e“ 


2t 


The last term is bounded as t j, 0, and the rest of the proof is as before. 


□ 
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